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The generation of non-gaussianity is studied in a three fluid curvaton model. By utilizing second 
order perturbation theory we derive general formulae for the large scale temperature fluctuation 
and non-gaussianity parameter, }nl, that includes the possibility of a non-adiabatic final state. In 
the adiabatic limit we recover previously known results. The results are applied to a three fluid 
curvaton model where the curvaton decays into radiation and matter. We find that the amount 
of non-gaussianity decreases as the final state of the system becomes more adiabatic and that the 
generated non-gaussianity in the scenario is small, |/jvl| ~ 0(1). 
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I. INTRODUCTION 

66 '. 

The Cosmic Microwave Background (CMB) has proven to be a feast of information for modern cosmology. Ever 
■ since its discovery by Penzias and Wilson in 1964 [l(, there have been a number of different experiments starting from 
RELIKT-1 0] resulting in drastic improvements to the quality of CMB data. A major scientific breakthrough was 
reached in 1992 when the COBE-satellite was able to detect the precense of the anisotropy in the CMB [!, [3, [H . 

The most recent data gathered by WMAP 0, Q is consistent with the hyphothesis that these perturbations were 
generated in an era of cosmic inflation. This is generally achieved with a slowly rolling scalar field which leads to an 
exponential expansion of the universe and the observed anisotropy is generated by the fluctuations of this inflaton 
field. Such a minimal scenario leads to adiabatic and gaussian perturbations in accordance with the current data. 
A well motivated alternative to the simplest inflationary scenario is the curvaton mechanism [T3, 
43 ; [II, 111 [II, [II S3, SIL SI SI S which the perturbations 

O |. are generated by a second scalar field dubbed the curvaton which stays subdominant during inflation and the actual 
JL ' expansion of space is still driven by the inflaton. This allows the inflation potential to have more natural properties 
compared to the single field scenario and can still lead to adiabatic perturbations. However, the extra degrees of 
freedom in the system now allow for the possibility that the final state is not necessarily purely adiabatic. Instead, 
the generation of an observable amount of isocurvature perturbations is a possibility that can distinguish the curvaton 
scenario from the simple single field inflationary model. 

Another ingredient that can differentiate the curvaton scenario from the simple inflaton hyphothesis is the concept 
of non-gaussianity, which has become more relevant with the data gathered by WMAP and with the advent of the 
' Planck satellite. The current limits from the WMAP 5-year data state that the local non-gaussianity parameter 
CO Jnl is limited to the values -9 < f^f 1 < 111 0. The single field inflationary scenario produces very little of 
non-gaussianity, Jnl ~ 0(1) whereas the curvaton scenario might lead to an observable non-gaussianity. 

The presence of non-gaussianity is therefore of paramount importance since it can differentiate between different 
CO 1 scenarios of the early universe. This paper is focused closely on the generation of non-gaussianity in a three-fluid 
model of curvaton decay where the curvaton decays into both radiation and matter. The concept of non-gaussianity 
in this scenario was first presented in [43| through first order perturbation theory. The study of non-gaussianities 
is however essentially dependent on the second-order terms and therefore the use of first order theory is not always 
justified. Our calculations incorporate the second-order perturbations from the start. Since the curvaton scenario 
might lead to adiabatic and isocurvature perturbations we have derived a generalized non-gaussianity parameter that 
includes the adiabatic state as a special case. 

This paper is organized as follows. In section II we present the relevant quantities and equations of motion of 
the perturbation up to second order. In this section we also derive the generalized non-gaussianity parameter /nl- 
In section III we present the equations of motion of the curvaton model and generalize the conserved quantity first 
presented in (43j. In section IV we present numerical results for /jvl in the three fluid curvaton model both for the 
case of constant and time-dependent interaction. We end this article with discussion and conclusions in section V. 
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II. PERTURBATIONS AT FIRST AND SECOND ORDER 



The theory of second order cosmological perturbations has been studied rigorously in the recent years. We will 
follow closely the notation of [§| and use a spatially flat Friedmann-Robertson- Walker -background. The metric tensor 
can be in this case expanded up to second order in the form Q 

g^daPda? = - (1 + 2^ + ^ 2 ))dt 2 + a(i)(^ (1) + ^-6j ( i 2) )dtdx i 

(1) 

dx l dx :> , 



where a(t) is the scale factor and 4^ r \ Ci}f\ and xlj are perturbation functions defined in 0] at first (r = 1) 
and second order [r = 2). Written in this form, different gauges can be straightforwardly given in terms of the 
perturbation functions: for example the Poisson gauge is defined as ui^ — — Xij — an d the spatially flat one 
is V (r) = x (r) =0 0. 

A useful set of equations can be derived from the continuity equations = Q", where Tf v is the energy- 
momentum tensor, Qi describes the energy transfer between different fluids and ; denotes the covariant derivative. 
From the continuity equations it follows that the equations determining the background evolution of individual fluids 
are 

Pi = -3H{l + co i )p i + Q i , (2) 

where u>i = Pi/ pi is the equation of state of th zth fluid and ' = d/dt i.e. derivative with respect to physical time. 
At first order one finds the evolution equations of the perturbed energy and pressure densities (on large scales) (4lj | 

5pf + 3H(6p^ } + - 3( Pl + P, ; )^ (1) = + SQ\ l) . (3) 

and at second order 

Spf + iH{5p {2) + «f } ) - 3( Pt + P,)^ (2) - H {1) [5 Pl + 5Pi + 2( Pi + P,)V (1) ] (4) 
= Qi^ + 5Q\ 2) Q^f + 2^6Q< i 1) . 

In addition to these, we also have the Einstein equations which can be used to give additional limits on the 
perturbation equations. For example in the Poisson gauge at first order ip^ = </>W and on large-scales 2ip^ 1 ' = —Sp/po, 
whereas in the spatially flat gauge tp^ = x^ = and 2cf)^ = —Sp/po @- 

At second order the equations get more complex and they are presented in detail in Q. The equations which we 
will need are the — and i — j-components of the Einstein equations on large scales which read in the Poisson gauge 
in a matter dominated universe as 

2 po 

^(2) _ (2) = _ 4^(13)2 _ i2 V - 3 (0ty (1) fiW (1) ) (5) 



In the spatially flat gauge the — 0-component is 

^=-^+4(0«) 2 . (6) 
2 pa 

A. Curvature perturbations 

An elegant way to study the evolution of perturbations is to use gauge- invariant curvature perturbations, which 
relate to curvature perturbations on homogeneous-density surfaces. At first order they are defined for component i as 

Cf ] = -<A (1) - (7) 



2 



where ' = d/d(ln(a)). At second order the corresponding quantity is defined as 

r (2) £ (!)/ : t 1 ) £ W r (!) 9 // 

P\ Da Da Da V Oa / \ 0- / 



Pi 



(8) 



Note that we are here neglecting gradient terms since we are only interested in the large scale behaviour of perturba- 
tions. 

The equation of motion of the first order curvature perturbations can be derived from eq. ([3]) and Einstein equations. 
The result is 



(1)/ _3&Pint(i) 



5Q 



int(i) 



where 5P int{t) = SP, -p'^pjp^ and SQ int ^ = 5Q % - Q'i^oJ X- 
At second order the corresponding equations read as [471 ] 
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(10) 



There are instances when the definition of different curvature perturbations might fail e.g. when p\ = 0. Therefore 
in our numerical evaluations we have used the spatially flat gauge and evaluated the density perturbations of different 
components in this gauge at first and second-order. The corresponding equations of motion can be easily read from 
equations Q and ([4} by going to the spatially flat gauge ip^ = x = 0. 



B. Non-gaussianity 



The generation of non-gaussianity in the two-fluid curvaton model at second order has been considered previously in 
in [jj [52j • In the two-fluid model the final state is adiabatic whereas in the three fluid model a significant isocruvature 
component is a possibility. Here we generalize the the results of Q to include the possibility of a non-adiabatic final 
state. We follow the standard notation presented in [4!| and use the Poisson gauge. Since we are interested only in 
the large scale non-gaussianities we can safely ignore some terms from the full expression J49(, including the integrated 
Sachs- Wolfe effect. The temperature fluctuations can in this approximation be written [9j as 



AT 



<,(2) 



(2)) _ 



l) r (l) 



Em 



(11) 



where <j)^ r > are the metric lapse functions, r 

= r (l) + 1 T (2) is 

the intrinsic fractional temperature fluctutation r 
AT/T\E m and all the terms are evaluated at the time of emission. 

We introduce new variables relating the final values of perturbations to their initial values: 



n 



/■(l)l 
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c (1) ' 
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VSfj.in/ 


92 " (C (1) ) 2 



(12) 



where the different numerators are evaluated at the time of decoupling, when the universe is matter dominated and 
£M ~ . The system is adiabatic if q\ = r% and q 2 = r 2 . 

(2) 

Since the universe is matter dominated during decoupling we can write Qn in the form 



c (2) = 

From the definitions of r\ and r 2 it follows that 

6 2) = r 2 (C {1) ? 



1 SpW 
3 Po 



5 f5p( 
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25r 2 
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(l)V 2 



(13) 



(14) 
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where we have also used the equation <jft-> = —3Q 1 '/5 which is valid on large-scales in a matter dominated universe. 
By combining these equations with the Einstein equations we can write (f>^ in a matter dominated universe in 
the form 



d2) 



16 
3 



5 r2 



(^ (1) ) 2 + 2v- 2 (av (1) 5^ (1) ) -6v- 4 (a i a J YflSi^ (1) a , > (i: 



(15) 



where the inverse Laplacians are to be understood as a formal expression. 

The intrinsic fractional temperature fluctuations t^ 1 "' can be written at first order as 



r (D = _ 



4 p 7 



X (!) 

Spy 
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En 



= 4^ 



3 ri 



(16) 



where we have used the definitions of Q , Cm i Qij r i an d 1 
is 



-3C'' 1 -' /5. At second order the corresponding variable 



15p 



4 Pl 



From the definition of Q we find 



ISp? 



4 p 7 



Em 



(17) 



(18) 



where we have also used the definitions of r% and g 2 . 

Combining all of the above expressions for perturbations and substituting into the equation for the temperature 
fluctuations (fTTj) one finally has 



AT 6ri — 5qi 



3r 3 



25(g 2 - qf) - 60 qi n + 96r 2 - 30r 2 



6ri — 5gi 



6ri(6ri — 5<?i 



3r x 



6ri - 5^i 



(19) 



where we have ignored the momentum dependent terms because we are interested in the large scale non-gaussianity 
@. To our knowledge this general formula has not been presented before. From eq. (|19p we can read that in the 
adiabatic limit, i.e. n = q\ and r 2 — q2, we recover at second order an extension of the first order Sachs- Wolfe effect 
AT/T = (j)^ /3 given in % 



AT 



-UP) - -(rhW) 



V" 



(20) 



In the opposite case of pure isocurvature perturbation, equation (|19p is not valid since we have assumed that r\ ^ 



and for isocurvature perturbations = 



Dec 



c, 



(1) 



Dec 



nC (1) 



We can now define the non-gaussianity parameter Jml in the general case. Following the notation of [51| we write 
the temperature fluctuations in the form 



AT 
~~f~ 



9t 



(i) 



NL\ 



cm 2 ' 



(21) 



where the factor gT depends on the state of the system e.g. for a completely adiabatic one gr — 1/3 and in our 
calculations gr = (6ri — 5qi)/(3r\). We have ignored the gradient terms from our definition of non-gaussianity but 
when calculating the bispectrum they need to be included. Note that this definition leads to a sign difference when 
compared to the usual approach using the Bardeen potential Q. Now from equation ()19p we can easily read the 
non-gaussianity parameter 



Jnl — 



25(q 2 



9?) 



60<?iri + 96r 2 - 30r 2 



6ri(6ri - 5qi) 



(22) 



which is well defined since r\ ^ and q\ <r%. 

In order to compare with observations we need to calculate the bispectrum, which requires that we take into account 
other effects e.g. the integrated Sachs- Wolfe effect and possible evolution after horizon crossing. This is beyond the 
scope of this paper and hence left for further work. For details of this procedure we refer to [H, [50| where the radiation 
transfer functions depend now on the adiabicity of the system via gT- 
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III. THE CURVATON MODEL 



A. Evolution equations 



We denote the curvaton by subscript (a) , radiation by (7) and matter by (to) . Since the curvaton field is oscillating 
it can be safely estimated [63j to behave like non-relativistic matter, i.e. oj a = 0. For radiation and matter we have 
uj 1 = 1/3 and ui m — and the interaction terms are [43| 



Qa = -T 1 f 1 (N)p a - Y m f m (N)p a 

<3 7 = T lfj( N )P<? 



(23) 



Qm = r m / m (iV)p CT , 



where I\ denote the strength of the interaction and functions fi(N) allow for time dependent interactions with 
iV = ln(a). 

Background equations @ can be written in terms of fractional densities Sli = pi/p for which the equations of 
motion are 1431: 



Hp 

Si' = Sl 7 {Sl 7 - 1) 4 



Q7 



(24) 



Hp' 



1 



3 



-ft 



From the definition of Sli it can be easily seen that Sl a + Sly + Sl m = 1, which means that one equation of motion of 
Sli is redundant. 

From eq. Q we can read the equations of motion for the first order density perturbations 
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in the spatially flat gauge. At second order the corresponding equations in the flat gauge are 

2 
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(26) 
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We also need the gauge invariant curvature perturbations Q in order to calculate the non-gaussianity parameter 
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Jnl- In the spatially flat gauge these are 



AX) = _ S P 



(i) 



-(2) _ 



Pi 

X ( 2 ) 
j 'Pi 

P'i 



2-3(1 + ^)] (Cf } ) 
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<5Q 
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Q- 1 Q l p' 



p[H 2 p^p 



(c 



(i) 



The set of equations P5|l and can now be evaluated numerically once the initial values have been set. 

We have chosen the system to be initially radiation dominated and non-adiabatic at first and second-order. The 
energy density of the curvaton field is during oscillations p a — m 2 a 2 where a is the amplitude of the field. A small 
perturbation of the field 5a leads to [321 ] 



Pa = Pa 

which means that at second order 



SpP + \tp { a ] + 0{5 P f) = m 2 (a% + 2a 5a + {5a? 



5p, 



(2) 



- ( 



(1) 



Pa <^ Pa 

This forces the second order perturbation of the curvaton to be 



a(2) = 



2 2H 7 3H 



2H(3H + T a ) 



(<?>y 



(28) 



(29) 



(30) 



which follows from eqs. (|29|) and (JSJ) after short calculations. Now since the first order curvature perturbations are 
linear equations we can scale them and set C^in = By pluggins this into eq. ([30]) and using the inequality T a <C Hq, 

1/2. The other perturbations are initally set to be zero 



-(2) 
^<7, in 



'•' ■ C7n = Cm'in = Cy"to = Ct'in = since we have assumed that only the curvaton field has initial perturbations. 



valid in our calculations, we can safely estimate Q 

_ /•(!) = /-(2) = >(2) 
^m, in ^7, in ^m.in 

Based on these initial values and equations ([24| , ([25| and ([26| we have calculated the amount of non-gaussianity in 
two different situations: (1) the curvaton decays into radiation and matter with constant couplings, i.e. f-y(N) = 1, 
f m (N) — 1, and (2) / 7 (iV) = 1 and the decay of the curvaton field into the matter component has explicit time 
dependence f m {N). We take f m (N) to be a continuous function with f' m {N) — at N = 0. 



B. Conserved quantities 

The authors of [43| noticed that the curvature perturbation ^iomp, which is related to a scaled matter like fluid 
component 

r 

Pcomp = Pm + -p — p Pa, (31) 

is a conserved quantity in the three-fluid model because 

r m 

Qcomp — Q in. H~ ~7~- j ~T~- Qcr (^^) 



and JQcomp = 0. This result can be generalised to all orders since <5Qcom P = and from equation (fTU|) we can easily 

(2) (2) 

see that this is equal to Qomp' = 0, i.e. Qomp is a conserved quantity. If we now include the possibility of time 
dependent interactions as in eq. (j2U|) . Qoinp is however no longer conserved because 

Q comp= ^ r ™y? ( ff m , p g ^0. (33) 



G 



Now in the both situations initially p r , 

Sp (D r M N) A „(i) 



(i) 



1 and Cmjn 



and we can see from equation l|27p that 



c 
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r m f m (N) § , 
r- r +r m f m (N)°' Ja 



-3(p„ 



r 7 +r m / OT (jv) 



r 7 +r m / m (iv) 



(i) 



JV=0 



o r m / m (jv) 



At >(1) 

3p CT ViE 



7V=0 



.A 

3,o 



N=0 



1, (34) 



where we have used the fact that T a f(N)/H 



T 7 +r m / m (jv)' 

■C 1 in our calculations. Once the curvaton has decayed completely, 



N=0 

Pa\Dec < /Ombec and Cm \dcc — Ccomp- Now if f m (N) = 1, ciomp are conserved and therefore (^J n = Cm |z?ec- In terms 
of r\ and r2 this means that r% = 1 and r2 = 1/2 which limits the possible values of the resulting non-gaussianity 
parameter (j2"2"]) significantly. 



IV. NUMERICAL RESULTS 
A. Continuous interactions, f-,(N) = 1, f m (N) = 1 
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FIG. 1: Evolution of the second order curvature perturbations in the curvaton-radiation-matter system with T 7 = 10 10 and 



T m = 10 , when (a) Q a0 
perturbation. 



10" 2 and (b) fl CT0 



10" 



Smaller initial curvaton density leads to a larger isocurvature 



The continuous case was first presented in [43| whereas in [4J] we studied the parameter space in more detail using 
physically motivated constraints. We also pointed out that the first order perturbation theory was not able to give 
limits on generated non-gaussianity. We now apply the second order theory results to calculate the amount of non- 
gaussianity. We have now CnP |d cc = C^Vn = 1 and Cm |d oc = C^in — 1/2- The evolution of C^ was explained in detail 
in [43| . If the curvaton fluid begins to dominate the system before it decays almost all of the radiation originates 



from the curvaton and C7 

>(2) 



Mi 



C 



(i) 



1. The same reasoning also applies at second order and during decoupling 



r (2)| 
1,7 I Dec 



Cf( n = 1/2 i.e. in terms of q\ and q%, q\ = \ and (72 = 1/2 and the system is adiabatic. If the curvaton 



decays before it begins to dominate the system, the curvature perturbations of the radiation fluid remain smaller and 

therefore lead to both adiabatic and isocurvature perturbations. This behaviour can be seen in the figures l(a)[|l(b) 

(2) 

where we have plotted the evolution of the second order perturbations Q i = er, 7,m. This also shows that the 

parameters q± and q2 are not independent i.e. a small absolute value of q\ means that \qi \ is also small whereas in 
the opposite case (72 = 1/2 when q\ = 1. 

In terms of the non-gaussianity parameter the above reasoning means that eq. (|22p can be written in the form 



NL 



25(g 2 -gi)-60gi+81 
6(6 - 5gi) 



(35) 



From this we can see that an adiabatic system, i.e. qi = 1, (72 = 1/2, gives Jml = 17/12 w 1.42. The maximum value 
Jnl = 2.73 corresponds to values qi = 0.438, 92 = 0.5 i.e. the system is non-adiabatic. 

For different interaction strengths the results are similar and in good agreement with the analytical approximation: 
the maximum values of non-gaussianity are close to the expected value and Jnl decreases as the system becomes 
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more adiabatic. In a previous paper we studied the same system in detail [4J| and found that the system becomes 
adiabatic as the interaction strengths decrease, hence weaker interaction strengths lead to less non-gaussianity. 

B. Time dependent interactions, / 7 (iV) = 1, fm(N) = 6(N - JV») 

The time dependent scenario was presented in detail in [44| , where we found only small differences compared to the 
continous case. Our choice for the interaction function f(N) is f(N) = (tanh((iV — N*)/t) + l)/2 where r = 10~5. 
Major alteration comes from the change of Hq — > ff* < Hq, which shifts different regions upward compared to the 
time independent scenario. We find a similar behaviour here. In terms of the curvature perturbation Ccomp it 1S no 
longer conserved because Q C om P ^ 0. However since the function f(N) reaches value 1 very quickly, f'(N) — is true 
almost everywhere and Ccomp is almost conserved. The value of the non-gaussianity produced in this scenario is thus 
very similar to the constant interaction case. 

V. DISCUSSION AND CONCLUSIONS 

In this paper we have studied the generation of non-gaussianity in the three fluid model of curvaton decay by means 
of second order perturbation theory. In the first part of this paper we introduced general formulae of perturbation 
theory, different gauge conditions and concluded it with a derivation of a general formula for the Sachs- Wolfe effect on 
large scales. It includes the possibility of a non-adiabatic final state and simplifies to the adiabatic formulas presented 
previously in . In the second part of this paper we applied this formula to the three fluid model of a curvaton decay 
and studied the generation of non-gaussianity in the temperature anisotropy of the CMB. We find that in general the 
amount of non-gaussianity produced in this scenario is small, Jnl ~ 0(1) both for constant and dynamical interaction 
strengths between the fluids. 

Our results indicate that the three fluid model leads to less non-gaussianity than the two fluid one. This is especially 
true if the curvaton field fails to dominate the system which leads to large non-gaussianity in the two fluid model 
because then Jnl — — 5/(4r) where r is small [a, H|. If the detected non-gaussianity is small one cannot, however, 
conclude that the three fluid model is responsible for it since such an observation is in agreement with the standard 
inflationary scenario /jvx = —1/2 j^, [37l |. One way to break this degeneracy is to compare the generation of gravity 
waves in the curvaton scenario [48[ with the standard inflation which usually produces much higher level of gravity 
waves. Another approach is to use the isocurvature: A non-adiabatic universe with a small amount of non-gaussianity 
would indicate the three fluid curvaton model as a natural candidate. 
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